Stability of multiquantum vortices in dilute Bose-Einstein condensates 
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Multiply quantized vortices in trapped Bose-Einstein condensates are studied using the Bogoli- 
ubov theory. Suitable combinations of a localized pinning potential and external rotation of the 
system are found to energetically stabilize, both locally and globally, vortices with multiple circu- 
lation quanta. We present a phase diagram for stable multiply quantized vortices in terms of the 
angular rotation frequency and the width of the pinning potential. We argue that multiquantum 
vortices could be experimentally created using these two expedients. 
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I. INTRODUCTION 

Bose-Einstein condensation in trapped atomic gases 
was realized and observed in 1995 |[ |j. Later, 
those pioneering experiments were followed by the cre- 
ation of singly quantized vortices Q and vortex lattices 
1^, H in such systems. Particularly after the creation 
of vortices in dilute atomic Bose-Einstein condensates 
(BECs), there has been pronounced interest in study- 
ing vortex configurations |^ due to their inherent con- 
nection with superfluidity. Especially, the stability of 
vortices in BECs has been subject to intensive research 
|,|,|l|,|l|,|l|,|ll,|l|,p,p. 

Singly quantized vortices in nonrotating, harmonically 
trapped BECs are locally energetically unstable within 
the Bogoliubov approximation in the sense that their 
excitation spectra contain anomalous negative-energy 
eigenmodes having a positive norm §, 0. If dissipa- 
tion in the system is not negligible, this local instabil- 
ity implies vortices to spiral out of the condensate [||. 
However, the anomalous modes are shifted to positive 
energies {w.r.t. the condensate energy) under sufficient 
rotation of the system 0. Furthermore, self-consistent 
finite-temperature theories predict singly quantized vor- 
tices to be locally stable even in the absence of external 
rotation In fact, a sufficient criterion for 

the existence of single-quantum vortices in dilute BECs 
seems to be their global stability, i.e., their having a lower 
free energy compared to the nonvortex state. This can 
also be quaranteed by adequate rotation of the system. 

Presently, many of the intriguing properties previously 
explored in conventional superfluid systems are being 
witnessed in dilute BECs. In addition to singly quantized 
vortices and vortex lattices, multiquantum vortices have 
been observed in thin films of superfluid "^He and in 
bulk rotating '^He-A |2^. In superconductors, they can 
be stabilized with the use of holes, antidots, or colum- 
nar defects as the pinning centers |2^, ^ . Contrary 
to the situation in Helium superfluids and superconduc- 
tors, multiquantum vortices remain unobserved in dilute 
BECs. This is presumably due to their generic instability 
against dissociation into an array of single quantum vor- 
tices [|l^, |3[ |2^. Consequently, an external rotation of 



the system would rather increase the number of vortices 
than nucleate a multiply quantized vortex. Nevertheless, 
in suitable anharmonic traps multiquantum vortices have 
been argued to be (globally) energetically favorable with 
respect to vortex lattices Q . 

In this paper, we demonstrate that by using vortex 
pinning ^] , one could stabilize vortices with multiple 
circulation quanta in harmonically trapped dilute BECs 
subject to an external rotation. We show that the mul- 
tiply quantized vortex states become locally stable and 
energetically advantageous when suitable external rota- 
tion and localized pinning potentials are employed. The 
required pinning of vortices could be realized by using 
an intense, blue detuned laser beam We argue that 
using this scheme, multiquantum vortices could be ex- 
perimentally created. 



II. MEAN-FIELD APPROXIMATION 

The low-temperature dynamics of the condensate is de- 
scribed by the time-dependent Gross-Pitaevskii equation 



ihdMr, t) = H(r)$(r, t) + g|$(r, t)\H{r, t). (1) 

The condensate wavefunction ^{r,t) is normalized ac- 
cording to / |<I>(r)pdr — N, where N is the total num- 
ber of particles. Above, the strength of interactions 
g = Anh^a/M as expressed in terms of the mass M of 
the atoms and the s-wave scattering length a. The effec- 
tive single-particle Hamiltonian 
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V2 + Kr(r) + Vpinir) + n-{rx ihV) (2) 



contains, in addition to the kinetic energy term, external 
trap and pinning potentials Vtr(r) and Vpin(r), respec- 
tively, and the angular momentum term, arising from 
the rotation of the system at the angular velocity fJ. 

Stationary solutions of the form <I'(r,t) — (/)(r)e~*''*/'' 
to Eq. (|l]) obey the time-independent Gross-Pitaevskii 
equation 



/.^(r) =7^(r)</)(r)+5|^(r)| 



(3) 
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where /i denotes the chemical potential determined by the 
normalization condition for the condensate wavefunction. 

The solutions of the stationary GP equation are mini- 
mizers of the energy functional 

E = y"[$*(r)H(r)$(r) + ||$(r)l4]dr. (4) 

By using Eq. (||), the energy of the condensate may be 
written in the form 

E = t^N-^-j \cj^[r)\''Av (5) 

in terms of the chemical potential and the condensate 
wavefunction. 

Adding a small perturbation 

(5$(r, t) = e[uq(r)e-'-^'*/'' + i;;(r)e*^''*/'']e-*^*/'' (6) 

to a stationary solution of Eq. (|^), reinserting it to the 
GP equation and neglecting terms of higher than linear 
order in e yields equations equivalent to the usual Bogoli- 
ubov equations ||29| 

Cuq{Y)+g(t>^{v)vq{Y) = EqUq{Y), (7a) 
Cvq{T)+g<l,*^{Y)uq{v) - ~E,v,{y) (7b) 

for quasiparticle amplitudes w^(r),u^(r), and eigenener- 
gies Eq\ here q labels the quasiparticle states. Above, 
C = 'H(r) — ^ + 2(7|0(r)p, and the quasiparticle 
amplitudes must obey the normalization /[[^^(r)^ — 

Ug'(r)p]dr — 6qq', manifesting the bosonic character of 
the excitations. The condition of local energetic stability 
of the solutions to the GP equation is that there exists 
no positive-norm quasiparticle excitations with negative 
energies Eg in the spectrum of Eqs. (|^). 

In what follows, we consider a Bose-Einstein conden- 
sate radially confined by a harmonic trapping potential 
Vtr(r) = jAfw^r^ in cylindrical coordinates r = {r,9,z). 
Here lu is the harmonic frequency of the trapping po- 
tential. The system is rotated in the plane perpendic- 
ular to the symmetry axis of the trap. We study the 
stability of rectilinear multiquantum vortex lines of the 
form (f>{r) — (j){r)e'^™^ , located along the rotation axis. 
The vorticity m determines the number of circulation 
quanta in the multiquantum vortices. We consider pin- 
ning potentials of the form Vpin(r) = yl0(i?pin — r), 
where denotes the unit step function and the amplitude 
A :i> ^ Q. For numerically solving Eqs. (^) and (0), we 
use computational methods similar to those described in 
Refs. @ [9|, 

III. ENERGETICS OF VORTICES 

In order to approach thermal equilibrium, a physi- 
cal system gravitates to the state that minimizes its 
free energy. Consequently, multiquantum vortices may 
be observed in experiments, provided that their free 



energies can be made lower than those of any other 
configurations — especially those of vortex arrays. How- 
ever, in the absence of pinning, the energy of a multi- 
quantum vortex is greater than that of a collection of 
singly quantized vortices with the same vorticity m. As 
a consequence, multiquantum vortices have not yet been 
seen in dilute, harmonically trapped BECs. 

In the case of a Bose-Einstein condensate without a ro- 
tating drive, the free energy minimizing state is vortex- 
free. However, when the condensate is rotated, its re- 
sponse is to acquire angular momentum by nucleating 
vortices. Pre-eminently, there exists a critical rotation 
frequency fie at which the energy of the condensate con- 
taining a singly quantized vortex becomes equal to that 
of a nonvortex state. At higher frequencies, vortices be- 
gin to nucleate from the edge of the condensate and to 
form arrays in the condensate . 

In the absence of the pinning potentials in the con- 
densate volume, there exists no global minima in the ef- 
fective potential felt by a vortex for rotation frequencies 
57 < fic [0- By the effective potential we mean the en- 
ergy of the system as a function of the position of the 
vortex in the condensate. Adding a localized pinning 
potential on the trap axis lowers locally the effective po- 
tential and can create a global minimum in the vicinity of 
the trap symmetry axis, thus enabling vortex occupation 
there. However, for < Sic, the condensate volume out- 
side the pinning potential would still remain energetically 
disadvantageous for vortices. The presence of circulation 
quanta in the system does not change the situation due 
to the repulsive interaction between vortices p3[ . 

In conclusion, for rotation frequencies f2 < fJc, states 
containing vortices in the region between the edge of the 
pinning potential and the condensate boundary are not 
globally energetically stable. In this paper we restrict the 
study to rotation frequencies < Sic, and thus do not 
consider the energetics of nonaxisymmetric vortex states. 
Extending the analysis consistently to frequencies SI > 
Sic would be an elaborate task due to the complications 
arising from different vortex array configurations, mutual 
interaction between vortices and spatial dependence of 
the vortex self-energy. 



IV. STABILITY OF MULTIQUANTUM 
VORTICES 

A. Global stability 

Mere rotation of the condensate is not sufficient for sta- 
bilizing multiquantum vortices in harmonically trapped 
Bose-Einstein condensates. In order to render multiquan- 
tum vortex structures globally stable, an additional pin- 
ning potential is required besides the rotating drive. 

The computed energies of axisymmetric vortex states 
containing m circulation quanta are shown in Fig. ^ for 
^pin = /im and i?pin = 6 /zm, as functions of the angu- 
lar rotation frequency SI. In the absense of the pinning 
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potential, the energy of a singly quantized vortex be- 
comes equal to that of a nonvortex state at the critical 
rotation frequency ilc . Rotation does not change the en- 
ergy of the nonvortex state m = because it has zero 
angular momentum. In the absense of the pinning po- 
tential, the multiquantum states are not the true energy 
minima for > flc, due to the dissociation instability. 
However, when a pinning potential is added, the stabi- 
lizing frequencies of multiquantum vortices are shifted 
below f2c (see the inset in Fig. where the nonaxisym- 
metric states are not energetically favored. 

Figure presents the energy differences AE{m) — 
E{m) — Eym — 1) between adjacent multiquantum states 
at Q = flc as functions of the radius of the pinning 
potential. The intersections AE{m) = in Fig. || de- 
fine the radii i?pin for which multiquantum configurations 
with m circulation quanta become globally energetically 
stable at 17 = f2c- In the inset, there are shown con- 
densate density profiles for vortices with m — 0, . . . ,5, 
stabilized by sufficiently wide pinning potentials. The 
rather large stabilizing values of i?pin compared to the 
core size of an unpinned vortex are partly due to the 
low rotation frequency applied, but also reflect conden- 
sate's rather strong tendency to distribute the vorticity 
uniformly throughout the system. 

The computed stability phase diagram for multiquan- 
tum vortices in terms of the angular rotation frequency 
and the radius of the pinning potential is depicted in 
Fig. ||. Each phase covers an area of the parameter space 
where a multiquantum vortex with given vorticity m is 
the minimum configuration of the energy functional. The 
phase boundaries are obtained by finding for each value of 
i?pin the rotation frequencies Q for which multiquantum 
vortices with successive vorticities have equal energy, see 
Fig. |l|. Specifically, the dashed fine in Fig. || corresponds 
to the vertical line at f2 = f2c in Fig. ||. 

Neither the rotation of the system nor the pinning po- 
tential alone suffices to stabilize multiply quantized vor- 
tices as is readily seen from Fig. || — a conclusion which 
applies equally well outside the parameter space covered. 
As the width of the pinning potential is increased, more 
circulation quanta 'fits' inside the pinning potential for 
a given rotation frequency. Similarly, as the rotation fre- 
quency is increased for a given pinning potential, the sys- 
tem can lower its energy by nucleating more circulation 
quanta in the multiquantum vortex located in the pinning 
potential. Evidently, the minimum radii of the pinning 
potential could be further decreased by increasing the an- 
gular rotation frequency above Qc — in other words each 
phase displayed extends beyond = fie, although the 
present study does not cover that regime. 



B. Local stability 

The remaining question is the local energetic stabil- 
ity of multiquantum vortices, which is determined by the 
sign of the lowest quasiparticle excitation energy of the 



condensate. For singly quantized vortices in nonrotating 
systems, there exists at least one anomalous mode, i.e., 
a negative-energy eigenmode with positive norm, within 
the Bogoliubov prescription. By rotating the system, 
those quasiparticle states may be lifted to positive ener- 
gies w.r.t. the condensate energy, implying local stability 
of the vortex state. Such stabilization by pure rotation is 
not possible, however, for unpinned multiquantum vor- 
tices as discussed below. 

The inset in Fig. ^ displays the energies of the low- 
est quasiparticle states for an unpinned doubly quantized 
vortex line. The system is rotated at the angular fre- 
quency fl/uj — 0.42, but the anomalous mode at qg — —2 
lies far below the condensate energy. Further increase of 
rotation results in a growing number of negative-energy 
excitations at higher values of angular momenta. The 
above analysis generalizes to higher values of m as well. 
And, the conclusion is that unpinned m > 1 vortices are 
locally unstable even in the presence of a rotating drive 

However, the use of a pinning potential changes the 
situation, and multiquantum vortices can indeed be also 
locally stabilized. In the main frame of Fig. ^ is shown 
the lowest excitation energies for stable m = 2 vortex at 
rotation frequencies JIl = 0.15 (o) and ily = 0.23 (•). 
The anomalous mode has disappeared from the spectrum 
as a consequence of the pinning, and the system is locally 
energetically stable even well beyond flu. The two-fold 
purpose of the pinning potential is thus to lift the anoma- 
lous vortex core modes to positive energies, and to lower 
free energies of multiquantum vortices, rendering them 
locally and globally stable, respectively. 



V. DISCUSSION 

In conclusion, we have studied multiply quantized vor- 
tices and their energetic stability in dilute Bose-Einstein 
condensates using the Bogoliubov approximation. Ener- 
gies of different multiquantum vortex configurations were 
computed and compared with each other in order to find 
the globally stable minimum-energy states. The analy- 
sis was restricted to rotationally symmetric states only 
by studying such rotation frequencies for which the nu- 
cleation of vortices outside the volume of the pinning 
potential is energetically hindered. We discussed both 
the local and global energetic stability of multiquantum 
vortex states and presented a phase diagram for their sta- 
bilization in terms of the radius of the pinning potential 
and the angular rotation frequency. 

In this paper, it was shown that a combination of a 
pinning potential and external rotation of the system fa- 
cilitates the existence of multiply quantized vortex states 
in harmonically trapped BECs. Such pinning of vortices, 
which could be accomplished, e.g., with an additional 
laser beam, has often been suggested but remains to be 
realized in the experiments |3|]. Based on the above 
analysis, we suggest that such multiquantum vortices 



could be created using available experimental techniques. Graduate School in Technical Physics for support. 

FIG. 1: Free energies per particle as functions of the angular 
frequency f2 for circulation quanta m = 0, 1, 2, 3, 4 and 5. The 
radii of the pinning potentials used are /um and 6 fim (in- 
set). In the absense of the pinning potential, the energy of a 
singly quantized vortex becomes equal to that of a nonvortex 
state at the critical rotation frequency flc/uj = 0.19. Notice 
the shifting of the intersections of the linos to lower values of 
n as -Rpin is increased. In the results presented in this paper, 
A^a/aho ~ 90, where aho = {h/ Mlu)^^^ is the harmonic oscil- 
lator length, a;/27r = 7.8 Hz, and N is the number of particles 
per harmonic oscillator length along the cylinder axis. 



FIG. 2: Energy differences l\E{m) = E(m) - E(m - 1) be- 
tween consecutive multiquantum vortex states at O = f2c 
as functions of the radius of the pinning potential -Rpin for 
m = 1,2,3,4 and 5. The points /\E(m) = determine the 
minimum radii of the pinning potential for which multiquan- 
tum vortices with m circulation quanta become energetically 
favorable. The subfigure shows stabilized, radial multiquan- 
tum density profiles for m = 0, 1, 2, 3, 4 and 5. The respective 
radii of the pinning potential Rpia{rn) are chosen to be (0, 0, 
6, 9, 12 and 15) /um. 



FIG. 3: Phase diagram for the stability of multiquantum vor- 
tices in terms of the angular rotation frequency Q, and the 
radius of the pinning potential i?pin. The parameter space 
is limited to frequencies < fie in order to guarantee the 
global stability of multiquantum states. Neither rotation of 
the system nor the pinning potential alone suffices to stabi- 
lize multiquantum vortices. Obviously, the minimum stabi- 
lizing radii of the pinning potential could be further lowered 
by increasing the rotation to frequencies O > Oc. However, 
the situation would be substantially more complicated in that 
regime due to the possible existence of stable states lacking 
rotational symmetry. 



FIG. 4: Lowest quasiparticle excitation energies for a stabi- 
lized {Rpin = 6 /jm) doubly quantized vortex line as functions 
of the quasiparticle angular momentum quantum number qe 
for rotation frequencies Ql = 0.15 (o) and ftu = 0.23 (•). 
The energies are measured relative to the condensate energy 
(■). For m > 1 the unpinned vortex carmot be made locally 
stable by simply rotating the system, as is demostrated by 
the inset for m = 2. The anomalous mode persists below the 
condensate energy and further increase of rotation frequency 
would bring about even more anomalous excitations for higher 
values of the angular momentum. 
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